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摘要 : A 2015 年 以 来 ， 人 们 已 经 探测 到 93 例 恒 星 级 双 黑 洞 并 合 引力 波 事件 ， 测 量 了 双 黑 洞 质 
量 、 自 旋 和 距离 等 参数 ， 并 获得 了 其 统计 分 布 。 不 同 双 黑 洞 起 源 机 制 产 生 的 双 黑 洞 性 质 分 布 ， 特 
别 是 自 旋 分 布 ， 会 显著 不 同 ， 其 可 用 来 限制 双 黑 洞 起 源 。 简 单 回 顾 了 双星 演化 和 动力 学 相互 作用 
等 双 黑 洞 形成 机 制 及 自 旋 来 源 ， 综 述 了 利用 自 旋 限制 双 黑 洞 起 源 的 一 般 方法 和 常见 模型 ， 以 及 利 
用 LIGO/Virgo 引力 波 数据 取得 的 相关 进展 ， 并 对 未 来 发 展 前 景 进行 了 展望 。 
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1 引 言 


LIGO (Laser Interferometer Gravitational-Wave Observatory) 于 2015 年 首次 探测 到 来 
自 于 宇宙 中 遥远 的 由 双 黑 洞 并 合 产 生 的 引力 波 事件 GW150914。 这 标志 着 人 类 不 仅 可 以 通 
过 电磁 信号 来 了 解 宇宙 ， 还 可 以 通过 引力 波 信 号 对 宇宙 进行 观测 ”“”。 引 力 波 的 观测 为 人 类 
认识 宇宙 提供 了 一 个 崭新 的 窗口 。 双 中 子 星 并 合 引力 波 事件 GW170817 及 其 电磁 对 应 体 的 
发 现 则 带领 人 类 走 进 多 信使 天 文学 的 全 新 纪元 。 从 此 ， 人 们 可 以 开始 结合 不 同 信使 从 不 同 角 
度 协同 研究 宇宙 中 重要 的 物理 过 程 ， 例 如 致密 星 的 形成 和 演化 、 中 子 星 的 内 部 结构 和 状态 方 
程 、 引 力 波 的 传播 等 ”。 

迄今 为 止 ， LIGO-Virgo 项 目 合作 组 已 经 在 激光 干涉 实验 中 探测 到 93 例 引 力 波 事件 ”。 

相 比 于 中 子 星 ， 黑 洞 的 质量 相对 较 大 ， 产 生 的 引力 波 信号 较 强 。 因 此 ， 绝 大 多 数 的 引力 波 
言 号 都 来 自 于 恒星 级 双 黑 油 (stellar-mass binary black holes， 的 并 合 过 
程 。 理 论 上 ， 这 些 双 黑洞 一 般 可 以 通过 以 下 几 种 机 制 产 生 : (1) 大 质量 孤立 双星 演化 机 制 
(evolution of massive binary stars， 即 EMBS channel); (2) 致密 zsh 力学 相互 作用 
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hannel); (4) 原初 黑洞 (primary black hole channel， 即 PBH channel)。 不 同 的 演化 机 制 产 
生 的 双 黑 洞 的 内 豪 物 理 参数 和 其 他 一 些 外 部 参数 分 布 不 同 ”“。 

引力 波 的 探测 为 我 们 提供 了 一 个 研究 这 些 机 制 的 独立 方法 。 双 黑洞 并 合 的 引力 波 时 序 
言 号 可 以 测 出 双 黑 洞 的 诸多 物理 参数 (包括 质量 、 自 旋 、 光 度 距离 等 )。 随 着 引力 波 探测 的 推 
进 ， 双 黑洞 并 合 样本 会 越 来 武大， 由 此 可 以 获得 双 黑 洞 的 物理 参数 和 一 些 外 部 参数 的 统计 分 
布 以 及 不 同 参数 间 的 关联 性 质 。 这 些 统计 性 质 将 会 极 大 地 帮助 我 们 理解 黑洞 这 种 宇宙 中 最 
奇异 的 致密 天 体 的 物理 性 质 和 双 黑 洞 的 形成 演化 。 通 过 分 析 最 新 的 03 引力 波 数据 ， 可 以 推 
导 较 准确 的 双 黑 洞 的 质量 和 自 旋 统 计 分 布 以 及 自 旋 - 质 量 关 系 等 ”“。 不 同 机 制 产生 的 双 黑 
洞 的 质量 和 自 旋 分 布 、 自 旋 - 质 量 关 系 等 可 能 存在 很 大 的 差异 。 目 前 ， 对 于 双 黑 洞 质 量 和 自 
旋 分 布 的 测定 和 分 析 将 可 以 对 其 起 源 机 制 进行 有 效 的 限制 ， 从 而 区 分 它们 并 限定 不 同 起 源 
对 观测 到 的 双 黑 洞 并 合 事 件 的 贡献 。 其 中 ， 利 用 双 黑 洞 自 旋 分 布 限制 双 黑 洞 起 源 成 为 近年 来 
的 热点 ， 已 有 一 些 工作 就 此 展开 了 深入 的 研究 “。 一 方面 ， 不 同 机 制 产生 的 双 黑洞 自 旋 分 
布 特别 是 等 效 自 旋 分 布 可 以 由 引力 波 探 测 样本 较为 准确 地 确定 ; 男 一 方面 ， 不 同 机 制 产生 的 
双 黑 洞 组 分 之 间 自 旋 的 相对 方向 分 布 差 异 明显 ， 理 解 上 简单 直观 。 

本 文 将 重点 对 恒星 级 双 黑 洞 的 自 旋 研 究 限制 其 起 源 机 制 进行 综述 ， 并 展望 未 来 发 展 。 本 
文 的 结构 如 下 : 在 第 2 章 将 简要 介绍 黑洞 的 基本 信息 ， 包 括 单 黑 洞 的 形成 过 程 以 及 黑洞 自 
旋 的 来 源 等 ;在 第 3 章 ， 我 们 将 介绍 恒星 级 双 黑 洞 的 不 同形 成 机 制 以 及 对 应 理论 预言 的 与 
自 旋 相关 的 信息 ; 第 4 章 介 绍 如何 通 过 自 旋 分 析 的 方法 限制 双 黑 洞 的 起 源 机 制 ， 包 括 贝 叶 
斯 参数 估计 方法 ，LIGO-Virgo 探测 到 的 致密 双星 并 合 样本 GWTC 数据 简介 ， 单 参数 和 多 
参数 的 自 旋 模型 等 ， 第 5 章 将 给 出 简要 的 总 结 ， 并 展望 未 来 双 黑 洞 自 旋 研 究 。 


机 制 (dynamical channel); (3) 活动 星系 核 / 大 质量 黑洞 协助 形成 机 制 (AGN/MBH-assisted 
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2 黑洞 基本 信息 介绍 


根据 广义 相对 论 ， 致 密 天 体 附 近 的 时 空 背 景 由 爱 因 斯 坦 场 方程 描述 : 


Guv = RT 5 (1) 


其 中 ，G 为 爱 因 斯 坦 张 量 ，T,, 为 描述 物质 分 布 的 能 动 张 量 ， 为 一 常数 。1916 年 ， 卡 
BLDG (Karl Schwazschild) 给 出 了 爱 因 斯 坦 场 方程 的 第 一 个 严格 真空 解 ， 即 史 瓦 西 度 
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2M 2M \~* 
ds? = (1 : Jae (1 ) dr? + r?(d6? + sin? 0d?) , (2) 


其 中 M 为 黑洞 质量 (G/c? = 1)。 史 瓦 西 度 规 在 + = 2M 处 有 奇异 性 ， 因 为 引力 红 移 效应 
的 存在 ， 观 测 者 永远 无 法 接收 到 这 一 位 置 的 光源 所 发 出 的 光子 。1958 年 ， 大 卫 . 芬 克 尔 司 坦 
(David Finkelstein) 将 > < 2M 的 区 域 解释 为 密度 如 此 之 高 ， 以 至 于 任何 东西 都 无 法 逃离 的 
区 域 。1967 F, INA EE (John A- Wheeler) 命名 该 区 域 为 黑洞 。 数 学 上 可 以 证 明 ， 克 
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尔 -纽曼 (Kerr-Newman) 黑洞 是 浙 近 平 直 稳 态 轴 对 称 黑洞 的 唯一 解 : 
A in? 0 
ds? = “2 (dt — asin? odo)? | 从 dr? + p*dé? 4 二 [@? Ha2)d0 adt]? ， (8) 


RH, =r? +a?cos0, A=r?—-2Mr +a +Q, a 为 黑洞 的 比 角 动量 。 惠 勒 等 人 进 一 
步 证 明了 黑洞 无 毛 定理 ， 一 个 渐 近 平 直 稳 态 黑洞 仅 由 黑洞 的 质量 M、 电 荷 Q、 自 转角 动量 
J = (aM) 3 个 参量 确定 。 在 天 体 物 理 当 中 ， 由 于 宏观 电 中 性 ， 我 们 通常 不 考虑 带电 的 黑洞 
模型 。 因 此 在 本 文中 ， 描 述 黑洞 的 最 主要 的 内 豪 物 理 参量 为 质量 和 自 旋 。 根 据 黑 洞 的 不 同 质 
量 ， 我 们 可 以 将 黑洞 分 为 : 恒星 级 质量 黑洞 、 中 等 质量 黑洞 和 超大 质量 黑洞 。 本 文 的 主要 研 
究 对 象 为 恒星 级 质量 黑洞 , 主要 通过 大 质量 恒星 演化 晚期 的 自 引力 黄 缩 形成 ， 同 时 可 能 会 伴 
有 超新星 爆发 和 /或 伽 马 暴 现象 。 

21 BE WAZ A 

宇宙 中 最 常见 的 天 体 是 处 在 不 同 演化 阶段 的 各 式 各 样 的 恒星 。 在 恒星 演化 晚期 ， 由 于 其 
核心 核反应 的 燃料 逐渐 消耗 殖 尽 ， 所 产生 的 辐射 压 不 足以 平衡 星体 自 引力 产生 的 向 内 的 压 
强 ， 恒 星 将 会 十 缩 形成 致密 星 ”。 对 于 质量 比较 小 的 恒星 (M < 9M。)， 将 会 形成 白矮星 。 
而 质量 稍 大 一 些 的 恒星 (Mo < M < 25M。) 在 其 演化 的 最 后 阶段 由 于 剩余 的 质量 比较 大 ， 
星体 的 自 引力 超过 电子 间 的 简 并 压力 ， 内 核 将 会 近乎 无 障碍 地 向 中 心 十 缩 。 而 博 缩 释放 的 引 
力 势 能 将 会 产生 极 高 的 温度 ， 从 而 导致 内 部 中 子 化 并 产生 大 量 的 中 微 子 。 如 果 超 新 星 爆炸 
是 非 球 对 称 的， 中 子 星 会 获得 踢 出 速度 ， 同时 中 子 间 的 简 并 压 将 抵抗 引力 坊 缩 。 在 这 种 情况 
下 ， 下 落 的 物质 冲击 中 子 星 表 面 ， 形 成 反弹 激 波 ， 最 终 将 中 子 星 以 外 的 外 壳 炸 开 ， 星 体 最 终 
将 会 演化 为 中 子 星 。 如 果 其 具有 一 定 的 自转 角 动 量 ， 则 会 形成 常见 的 脉冲 星系 统 。 

而 对 于 质量 更 大 的 恒星 ， 核 心 的 中 子 简 并 压 仍 然 不 足以 支撑 自 引力 的 拥 缩 ， 整 个 星体 
被 进一步 “无 限制 ”压缩 ， 最 终 形成 天 体 物理 学 常见 的 黑洞 。 此 类 黑洞 的 质量 一 般 大 于 中 
子 星 的 最 大 质量 ， 也 即 奥 本 海 默 静态 极限 (通常 认为 是 2.4M。， 但 实际 上 静态 假设 通常 不 
成 立 ， 仅 能 给 出 (2 ~ 3)Mo 的 范围 )。 由 于 特别 大 质量 的 前 身 星 自身 过 大 的 星 风 损失 以 及 演 
化 晚期 可 能 存在 的 正 负 电子 对 不 稳定 性 (pair instability) 等 原因 ， 人 恒星 级 黑洞 的 质量 通常 被 
限制 在 50/60Me 以 下 。 值 得 注意 的 是 ， 根 据 恒星 演化 理论 的 模拟 结果 显示 ， 人 恒星 级 黑洞 在 
(50 ~ 120)Mo 附近 存在 质量 空缺 ( 见 图 四 。 通 过 电磁 波 观 测 ， 人 们 发 现 并 给 出 了 银河 系 中 
二 十 多 颗 黑 洞 的 质量 ， 大 多 数 分 布 在 (5 ~ 30)M 范围 内 "“。 然 而 ， 通 过 引力 波 的 双 黑 洞 并 
合 观 测 发 现 了 一 些 双 黑洞 的 主 黑洞 质量 可 以 大 于 50/60Ms， 对 现 有 的 黑洞 形成 模型 提出 了 
Pie” 

2.2 ” 单 黑洞 的 自 旋 

黑洞 的 自 旋 通常 来 自 于 前 身 恒星 的 角 动量 。 根 据 恒 星 不 同 的 角 动 量 耗 散 机 制 ， 将 会 产 
L 有 不 同 自 旋 的 黑洞 。 一 种 可 能 的 耗 散 机 制 是 通过 经 线 方向 的 物质 流动 以 及 满 流 的 相互 
用 来 转移 恒星 内 部 的 角 动量 。 如 果 这 个 过 程 中 有 磁场 参与 ， 根 据 Tayler-Spruit 磁 发 电机 
(magnetic dynamo) 机 制 ， 恒 星 角 动量 的 耗 散 将 会 变 得 更 加 有 效 。 
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金属 丰 度 (大 致 对 数 刻 度 ) 


9 10 25 34 40 60 100 140 260 
初始 质量 / M。 


图 1 ”大 质量 孤立 恒星 死亡 形成 致密 星 类 型 与 恒星 初始 质量 和 初始 金属 丰 度 的 关系 


在 引力 波 被 探测 到 之 前 ， 人 们 通常 利用 K 射线 波段 的 观测 对 双星 系统 中 的 黑洞 无 量 纲 
自 旋 a 进行 测量 。 在 双星 系统 中 ， 伴 星 的 物质 将 会 被 黑洞 吸 积 从 而 形成 吸 积 盘 。 通 过 对 吸 
积 盘 进 行 观测 ， 将 可 能 获得 黑洞 的 相关 参数 。 根 据 广义 相对 论 和 吸 积 盘 理 论 的 预言 ， 黑 油 
的 吸 积 盘 内 半径 会 被 截断 至 最 内 稳定 圆 轨道 处 (innermost stable circular orbit, ISCO), mi 
这 个 轨道 的 半径 Risco 是 黑洞 无 量 纲 自 旋 参 数 a = J/M? 的 单调 函数 ( 见 图 团 。 因 此 通过 
拟 合 吸 积 盘 的 连续 谱 还 原 吸 积 盘 的 尺寸 ， 可 以 测量 黑洞 的 自 旋 ”"。 而 另外 一 种 常见 的 方法 ， 
则 是 利用 Fe Koc 发 射线 由 于 引力 红 移 等 效应 产生 的 展 宽 估 计 吸 积 盘 内 边缘 的 位 置 ， 从 而 测 
量 黑洞 自 旋 “ 所。 除 此 之 外 ， 利 用 吸 积 盘 的 X 射线 强度 的 准 周期 振荡 也 可 以 获得 ISCO 相 
关 的 信息 ， 从 而 测量 黑洞 的 自 旋 。 但 这 种 方法 对 震荡 模式 的 模型 依赖 性 很 强 ”。 

以 上 几 种 方法 都 有 着 各 自 的 问题 。 例 如 ， 拟 合 连续 谱 需 要 对 吸 积 盘 内 边缘 的 倾角 进行 
测定 ， 往 往 其 误差 特别 大 ，Fe Ka 发 射线 展 宽 法 则 要 求 对 产生 硬 X 射线 的 电子 晕 进 行 精确 
建 模 ， 准 周期 振 莲 方法 要 求 对 吸 积 盘 的 振荡 模式 有 较为 精准 的 研究 ， 等 等 。 此 外 ， 迄 今 为 目 
能 够 使 用 这 些 方法 来 进行 观测 的 仅 20 余 个 源 ” ， 相 对 于 引力 波 观 测 而 言 并 不 多 。 因 此 ， 通 
过 X 射线 双星 对 黑洞 自 旋 进 行 测量 存在 局 限 性 。 通 过 对 引力 波 观测 到 的 大 量 双 黑 洞 并 合 事 
件 的 参数 进行 拟 合 ， 我 们 可 以 较 精确 地 测量 黑洞 的 自 旋 。 
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2 ”最 小 稳定 圆 轨道 半径 Risco 和 视界 面 半径 Ru 与 黑洞 自 旋 参 数 a. 的 关系 


3 恒星 级 双 黑 洞 形成 机 制 


如 果 两 个 恒星 级 黑洞 沿 着 准 开 普 勒 轨道 相互 绕 转 ， 将 会 形成 引力 自 束缚 的 恒星 级 双 黑 
洞 系 统 。 在 旋 进 的 过 程 中 ， 由 于 不 断 地 辐射 引力 波 ， 它 们 的 轨道 能 量 将 会 逐渐 降低 ， 因 此 两 
个 黑洞 逐渐 靠近 ， 并 最 终 经 历 并 合 以 及 后 续 的 铃 宕 阶段 ， 形 成 并 合 后 的 单 黑洞 。 例 如 入 类 
所 发 现 的 第 一 个 引力 波 事件 GW150914， 是 由 质量 分 别 为 35.61Mo 5 30.6Mo 的 恒星 级 黑洞 
并 合 产生 。 而 一 个 特别 的 引力 波 事件 GW190521 在 99.68% 的 置信 和 度 下 ， 两 个 黑洞 的 质量 
都 大 于 65Mo” ， 处 于 恒星 演化 数值 模拟 预言 的 质量 间隙 之 中 。 这 可 能 意味 着 :这些 质 量 在 
(50 ~ 120)Mo 的 黑洞 主要 由 初代 黑洞 并 合 产生 ， 而 不 是 直接 由 大 质量 恒星 演化 而 来 的 。 
理论 上 ， 目 前 主要 有 以 下 两 种 机 制 可 以 产生 恒星 级 双 黑 洞 : (1) 大 质量 孤立 双星 演 
化 机 制 …”，(2) 动力 学 相互 作用 机 制 (dynamical interactions in dense stellar clusters, 
Dynamical channel) 2029) | 不同 机 制 产 生 的 双 黑 洞 的 质量 、 自 旋 ， 以 及 轨道 特性 等 在 统计 学 
意义 上 各 有 其 特点 。 本 章 将 着 重 于 介绍 这 两 种 最 主要 的 双 黑 洞 形成 机 制 并 对 其 进行 讨论 ， 而 
对 其 他 可 能 产生 双 黑 洞 的 机 制 ， 例 如 活动 星系 核 以 及 大 质量 黑洞 机 制 (AGN/MBH assisted 
channel)” 和 原初 黑洞 机 制 ”” 仅 做 简单 的 介绍 。 
3.1 ”大 质量 孤立 双星 演化 机 制 

大 质量 孤立 双星 的 演化 机 制 通常 指 的 是 孤立 (不 与 环境 发 生 相互 作用 ) 双星 经 过 公共 包 
层 阶 段 最 终 演化 形成 双 黑 洞 的 过 程 。 该 理论 认为 ， 当 主星 率先 结束 其 主 序 阶段 的 HH 燃烧 ， 它 
的 He 核 将 会 发 生 收缩 ， 而 外 面 的 H 包 层 则 开始 膨胀 。 当 洛 希 办 被 完全 占 满 时 ， 来 自 于 伴 
星 的 潮汐 引力 将 会 对 抗 主星 的 自 引力 ， 从 而 使 得 物质 从 主星 转移 到 伴星 上 ”。 最 终 ， 主 星 
将 会 形成 一 颗 裸 露 的 He 核 星 ， 并 通过 星 风 不 断 损失 物质 ， 直 到 当 核 能 消耗 殖 尽 时 ， 主 星 将 
会 直接 夫 缩 为 黑洞 。 如 果 出 生 踢 出 (natal kick)?” 或 者 物质 损失 不 能 破坏 这 样 的 黑洞 -伴星 系 
统 ， 那 么 他 们 将 进一步 演化 ， 直 到 伴星 开始 膨胀 并 且 其 物质 被 吸 积 到 黑洞 上 。 由 于 最 开始 伴 


Gi 
Gi 
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星 已 经 从 主星 吸 积 了 一 些 物质 ， 以 及 主星 本 身 的 物质 损失 ， 可 能 的 星 风 ， 超 新 星 爆发 等 ， 在 
这 个 阶段 伴星 的 质量 将 会 远 远大 于 主星 的 质量 。 该 系统 的 物质 转移 过 程 将 变 得 十 分 不 稳定 ， 
这 种 动力 学 非 稳定 性 使 得 双星 进入 公共 包 层 演化 阶段 。 当 黑洞 和 He 核 组 成 的 内 双星 与 公共 
包 层 之 间 存 在 较 差 转动 时 ， 拖 电 的 力量 将 会 使 得 双星 轨道 能 量 发 生 耗 散 ， 从 而 可 能 引起 公共 
包 层 的 物质 发 生 抛 射 。 在 这 个 阶段 ， 伴 星 将 会 失去 自己 的 包 层 变 成 另 一 颗 裸露 的 He 星 ， 最 
终 演化 圾 缩 形成 一 个 黑洞 。 如 果 假 设 第 二 次 超新星 爆发 没有 破坏 这 个 双 黑 洞 系统 ， 它 将 会 进 
一 步 经 历 旋 进 ， 并 合 ， 以 及 铃 宕 阶段 ， 最 终 产 生 新 的 较 大 质量 的 黑洞 ”( 见 图 团 。 这 其 中 每 
一 个 过 程 产生 的 引力 波 信号 都 可 能 被 引力 波 探测 器 所 探测 到 。 需 要 指出 的 是 ， 以 上 机 制 形成 
双 黑 洞 的 过 程 中 ， 物 质 转移 和 公共 包 层 演化 等 物理 过 程 非常 复杂 ， 目 前 对 其 理解 有 很 大 的 不 
确定 性 ， 因 而 可 能 会 影响 结 


t/Ma a/Ro e 


0.0000 | MS 96.2Mo 60.2Mo |2463 0.15 


3.5445 | HG 92.2Mo 59.9Mo |2140 0.00 


3.5448 | 或 42.3Mo 84.9Mo |3112 0.00 


3.8354 | He 星 ”39.0Mo 84.7Mo |3579 0.00 


3.8354 | BH 35.1Mo 84.7Mo | 3700 0.03 


5.0445 |BH 35.1Mo 82.2Mo |3780 0.03 


5.0445 |BH = 36.5Mo & 36.8Mo | 43.8 0.00 
5.3483 |BH 36.5Mo è Ba He 星 34.2Mo| 45.3 0.00 
县 直接 黑洞 
5.3483 |BH 36.5Mo Pe BH 30.8Moe | 47.8 0.05 
县 合并 
10.294 © 0 0.00 
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由 于 双星 各 子 星 的 自 旋 较 快 ， 动 力学 的 不 稳定 混合 造成 了 化 学 类 均匀 演化 。 这 种 演化 过 
程 中 ， 恒 星 几乎 不 会 膨胀 ， 星 风 也 很 小 ， 形 成 的 He 星 质量 很 大 ， 因 此 可 以 形成 较 大 质量 的 
黑洞 ， 具 体 解 释 请 查阅 文献 [34 中 简介 部 分 第 3 段 。 图 田 简 要 说 明了 化 学 类 均匀 演化 与 一 
般 恒星 演化 的 区 别 。 


`‘ 
Se = 
3a A” TDN 3-b 洛 希 办 溢出 
4 *% 
I \ 
| @ > 
\ a 
\ / 
Se we 
bea eo 并 合 或 通过 洛 希 准 有 滋 出 
造成 严重 质量 损失 
| ) 
\ pa 
2 
NSS we 


图 4 (828 5A aE eR 


另外 ， 极 端 贫 金属 的 第 一 代 恒 星 ( 星 族 TT) 也 可 以 产生 双 黑 洞 。 在 此 情况 下 ， 由 于 前 身 
双星 中 一 开始 不 含 C、N、0O 等 元 素 ，H 聚变 只 能 通过 质子 -质子 过 程 产生 。 当 恒星 核 温 度 
升 高 至 可 以 发 生 3-a 反应 之 后 ，CNO 循环 才 开 始 主导 核心 的 核反应 。 这 样 的 星 族 III 恒星 
演化 后 期 壳 层 的 膨胀 被 抑制 ， 从 而 避免 了 进行 物质 交换 。 最 后 ， 两 颗 恒 星 各 自 替 缩 ， 从 而 形 
RAE. 

简单 的 理论 计算 给 出 了 上 述 公 共 包 层 演化 ， 化 学 均匀 演化 以 及 星 族 III 演化 模型 在 本 
地 > = 0 的 并 合 率 ， 分 别 约 为 50、10、0.02 Gpec 习 .a-1 ”9 ， 而 这 些 黑洞 的 典型 质量 为 
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(20 ~ 80) Mo 


在 这 些 不 同 的 大 质量 孤立 双星 的 演化 模型 当中 


的 自 旋 与 双星 的 轨道 角 动 量 平行 。 如 


FP， 一 般 都 假设 双星 系统 中 两 颖 前 身 恒星 
黑洞 不 经 历 出 生 足 出， 那么 单纯 的 物质 交换 和 测 汐 相 


41 卷 


互 作用 将 不 会 显著 地 改变 它们 的 自 旋 方 向 。 这 样 直接 导致 了 形成 的 双 黑 洞 系统 两 个 黑洞 的 


自 旋 倾向 于 平行 。 另 一 方面 ， 如 果 黑 洞 的 质 


=. 
E 


出 生 踢 出 ， 自 旋 方向 将 会 变 得 不 平行 。 


3.2 ”动力 学 相互 作用 机 制 


动力 学 相互 作用 机 制 是 
参与 引力 作用 都 可 以 被 归结 到 这 种 产生 机 人 


上 中。 一 般 地 ， 我 们 可 


比较 小 ( 约 小 于 10Mo) 且 存 在 超新星 爆发 的 


的 说 法 ， 只 要 在 形成 双 黑 洞 的 过 程 中 有 其 他 星体 
以 根据 动力 学 过 程 发 生 的 环 


境 进 行 细 分 ， 例 如 球状 星团 和 核 星 团 中 的 演化 ， 等 级 三 体系 统 中 的 演化 等 。 


在 球状 星团 中 ， 初 始 单 黑洞 的 质量 相对 较 大 ， 将 会 | 


域 。 一 旦 到 了 如 此 密集 的 区 域 ， 黑洞 将 可 能 立即 形成 
过 程 参 与 到 其 他 的 黑洞 -恒星 系统 
核 星 团 中 ， 一 般 认为 存在 一 个 超大 质量 的 
质量 分 层 效 应 更 加 明显 ， 将 提供 更 有 利于 双 黑 洞 产生 

另外 一 个 常见 的 动力 学 机 制 便 是 宇宙 中 
内 双星 和 较 外 的 第 三 颗 黑 洞 将 会 发 生 角 动 
椭 率 被 激发 到 非常 高 的 值 时 ， 内 双星 的 轨 ; 
促进 他 们 发 生 并 合 Kozai-Lidov 效应 。 在 密 
的 双 黑 洞 自 旋 通常 被 认为 是 方 


将 较 轻 的 
F ( 见 图 回 。 


于 质量 分 层 效 应 掉 入 星团 的 中 心 区 


自 束缚 的 双 黑 洞 系统 ， 抑 或 是 通过 三 体 
体 踢 出 ， 形 成 新 的 双 黑 洞 系统 。 而 在 


它 将 会 产生 巨大 的 引力 势 场 ， 使 得 


向 随机 的 。 这 是 因 


的 高 密度 环境 。 

上 分 常见 的 等 级 三 体系 统 ”。 在 这 样 的 系统 中 ， 
转移 ， 从 而 使 得 内 双星 的 椭 率 不 断 变化 。 当 
各 可 以 通过 引力 波 或 者 测 汐 力 耗 散 掉 ， 从 而 
场 或 是 三 体系 统 中 通过 动力 学 相互 作用 产生 
此 复杂 和 随机 的 动力 学 相互 作用 过 


程 中 ， 我 们 很 难 在 物理 意义 上 找到 一 个 特别 的 方向 作为 黑洞 自 旋 的 初始 方向 ， 所 以 一 般 都 认 
为 由 此 机 制 形成 的 并 合 双 黑洞 系统 初始 自 旋 是 完全 随机 的 ， 而 根据 后 牛顿 轨道 的 数值 计算 


显示 ， 这 样 的 随机 性 将 会 被 保持 到 双 黑 洞 系统 最 终 } 
制 不 同 ， 动 力学 机 制 产生 的 双 黑 洞 自 旋 方 向 通常 认为 完全 随机 。 


3.3 ”其 他 可 能 的 机 制 


除了 以 上 两 种 最 主要 的 形成 机 和 
洞 协助 机 制 、 原 初 黑洞 机 外 


Ifo 
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定性 直接 产生 的 恒星 ， 并 通过 吸 积 成 长 为 更 大 质量 的 恒 
的 相互 作用 ， 形 成 大 质量 双星 和 双 黑 洞 系统 ”。 而 原初 黑 》 


=I 


EE 级 质量 双 黑 洞 也 可 能 
a 


此 ， 与 大 质量 孤立 双星 演化 机 


通过 活动 星系 核 / 大 质量 黑 


P 心 黑洞 周围 存在 一 个 巨大 的 吸 积 盘 。 


i 获 一 些 恒星 ， 抑 或 是 通过 盘 不 稳 
这 些 恒星 通过 与 周围 气体 和 恒星 


= 


有 学 者 认为 相当 一 部 分 上 暗物质 是 由 


以 产生 恒星 级 双 黑 洞 的 并 合 引 力 波 信和 号 


司 来 自 于 早期 宇宙 的 密度 涨 落 ”。 
级 质量 的 原初 黑洞 构成 ， 它 们 之 间 发 生 的 并 合 也 可 
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图 5 
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抛射 时 间 :2.7Ga T ane AU 
a=0.38AU 到 旋 进 阶段 时 间 :8.7 Ga a 
e=0.76 


抛射 时 间 : 4.2 Ga 
到 旋 进 阶段 时 间 : 6.3 Ga 
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4 基于 目 旋 分 析 对 双 黑 洞 起 源 的 限制 


ppl 


a, LAMA 


其 
对 
万 
叶 


=I 


自 旋 的 相关 信息 将 会 提供 关于 恒星 级 双 黑 洞 形成 机 制 的 重要 信息 。 
洞 的 有 效 自 旋 Xeg 十 分 敏感 ，Xeg ENF”: 
1 


Xeff = M 
HP, my m 分 别 为 双 黑 洞 的 质量 ，A4 = mi 十 m2 为 质量 
角 动 量 (假设 为 z FTI), si. s2 为 双 黑 洞 相 应 的 无 量 纲 黑洞 自 旋 ， 
cos 61. cos 62 是 自 旋 角 动量 与 轨道 角 动 量 上 的 夹 角 。 在 旋 进 阶段 ， 有 效 
同 的 进 动 效 应 ， 引 入 一 个 新 的 参数 ， 有 效 进 动 自 旋 ， 定 义 如 下 : 


为 了 进一步 评估 双 黑 ; 


洞 的 轨道 


LI 


Mı mg 


a Al, Si, 


L 1 
(=: | 2) . = M (m1sicos01 + m2s2C0s02) , 


41% 


而 引力 波 波形 对 双 


(4) 


So 为 黑洞 的 自 旋 角 动 


中 ， = 2+4 3q/2 Seon q AMS) 


Xp = Bim? ee (Bi S11, Bo S| L) > 0 ; 


同 的 质量 比 。 我 们 可 以 利 月 


自 旋 几 乎 是 守恒 的 。 


(5) 


引力 波 的 探测 


双 黑 洞 与 上 述 两 个 自 旋 相关 的 参数 进行 测量 ， 从 而 基于 自 旋 分 析 对 恒星 级 双 黑 洞 的 起 源 ， 
本 章 将 主要 介绍 利用 贝 


其 是 对 大 质量 孤立 双星 演化 机 制 和 动力 学 机 制 进行 一 定 的 限制 。 
斯 统计 处 理 GWTC 引力 波 观 测 数据 的 一 般 方法 ， 
星 级 双 黑 洞 形成 机 制 的 研究 进展 。 


4.1 GWTC 数据 简介 
GWTC (Gravitational Wave Transient Catalog) 是 LIGO/Virgo/KAGRA 项 目 合 作 组 


公布 的 引力 波 探测 数据 集 。 


阶 


(O1) 


个 


双 


以 


者 


Xer 间 有 具有 很 强 的 关联 
(2) GWTC-2 (包括 
期 间 (2019 年 4 月 1 日 
体 。 相 比 GWTC-1 日 


4 


行 
选 


全 部 为 双 黑 洞 并 合 事件 


[40， 多， 


(1) GWTC-1 的 数据 包 全 LIGO 和 Virgo 自 2015 年 9 H 12 


[10] 


, 2016 £11 H 


黑洞 的 总 质量 范围 为 18. 
图 加 展示 了 推断 出 的 引力 波源 的 有 效 自 旋 Xefg， 以 及 有 效 进 动 自 旋 xp 的 后 验 分 布 。 可 
都 在 0 左右 ， 仅 有 GW170729 和 GW151226 两 个 源 显 
而 进 动 自 旋 x 


30 


看 出 几乎 所 有 双 黑 洞 的 Xeg 的 峰值 
ij 离 0， 分 别 为 : 0.11 ~ 0.58 和 0.06 ~ 0.38 (90% 置信 区 间 )。 
性 。 


口 


日 至 2017 年 8 H (02) 的 观测 结果 。 主 要 探测 至 


个 双 中 子 星 并 合 事件 、GW150914 以 及 GW151012 等 10 AI yF 


以 及 相应 的 单 参 数 和 多 参数 模型 限制 恒 


前 进行 了 三 次 发 布 ， 对 应 探测 器 运行 的 Ol 和 O02、0O3a、0O3b 


段 ， 共 有 93 个 确证 的 引力 波 事件 ， 其 中 有 两 例 中 子 星 并 合 ， 一 例 黑 洞 -中 子 星 并 合 ， 剩 下 


#4 2016 年 1 月 19 H 
包含 GW170817 一 
合 事件 。 其 中 探测 到 的 


6132Mo ~ 84.4+193Mo， 距 离 范 围 为 320+139 ~ 2840414 Mpc, 


的 数据 ， 双 黑洞 的 总 质 
图 四 展示 了 GWTC-2 探测 到 的 引力 波源 的 物理 参数 以 及 不 确 
共有 非 零 的 xeg， 且 在 90% 的 置信 度 情 况 下 ，Xeg 均 大 于 0。 其 中 具有 最 大 有 效 自 旋 的 双 黑 


2019 年 10 月 1 日 ) 的 观测 数据 ， 探 测 
量 范围 扩大 为 (14 ~ 182) Mo。 


» 与 有 效 自 旋 


2.1) 的 数据 包括 Advanced LIGO 以 及 Advanced Virgo Æ O3a 运 


到 44 个 新 的 引力 波源 候 


定 度 ， 其 中 有 


?这 里 采用 与 相关 文献 中 一 致 的 


旋 标 识 ，S 与 前 文中 的 J 均 代 表 自 旋 角 动量 。 


13 个 系统 
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YE: a) GWTC-1 结果 中 的 有 效 自 旋 参 数 后 验 分布 ，b) 白色 部 分 代表 GWTC-1 中 有 效 进 动 自 旋 参 数 的 先 验 分 
布 ， 彩 色 部 分 表示 GWTC-1 中 有 效 进 动 自 旋 参数 的 后 验 分 布 。 


[40] 


图 6 GWTC-1 观测 到 的 自 旋 参 数 数据 分 布 


IRL GW190517-055101, xer = 0.52+9:19。 进 动 自 旋 同 样 显 示 出 与 有 效 自 旋 xer 具有 很 强 
的 关联 性 。 
(3) GWTC-3 的 数据 包括 Advanced LIGO 以 及 Advanced Virgo 在 O3b 运行 期 间 (2019 
年 10 月 1 日 一 2020 年 3 月 27 A) 的 观测 数据 ， 增 加 了 O3b 期 间 探测 到 的 38 个 新 的 引力 
波源 或 候选 体 (包含 第 一 例 黑洞 - 中 子 星 并 合 事件 )， 参 数 范围 与 GWTC-2 的 结果 类 似 ， 见 
B 大 多 数 引 力 波源 的 有 效 自 旋 都 接近 于 0， 有 4 个 大 于 0 (置信 度 大 于 89%)。 但 值得 注 
意 的 是 ，GW191109-010717 和 GW200225-060421 在 90%, 85% 置信 区 间 里 ， 有 效 自 旋 小 
于 0。 进 动 自 旋 与 有 效 自 旋 Xeg 关联 性 变 得 更 强 。 
4.2 ”基于 贝 叶 斯 理论 的 分 析 方 法 简介 
在 引力 波 数 据 分 析 中 ， 由 于 引力 波 的 强度 h(f) 本 身 非常 低 ， 而 噪声 却 很 高 ，h(f) 的 范 
Fy 10-22 ~ 10-??， 我 们 通常 利用 模板 匹配 的 方法 对 引力 波 信号 及 噪声 进行 分 析 ”。 其 信 
品 比 可 以 通过 如 下 表达 式 计算 : 
fmax 2 
p=4| RO ap ， (6) 


min nf) 


cH, A) 为 引力 波 探测 器 的 单 侧 功 率 谱 ，fiwis 和 fax 为 引力 波 频率 的 积分 上 下 限 。 通 
常 我 们 取 p = 8 为 引力 波 探测 的 阔 值 (等 价 于 10% 的 误 报 率 )， 也 即 高 于 该 值 就 定义 为 探测 
到 一 个 引力 波 事件 。 

一 旦 引力 波 信号 被 接收 到 ， 那 么 将 可 以 利用 贝 叶 斯 统计 的 办 法 对 引力 波源 物理 参数 O 
的 后 验 分 布 P(O|M) 进行 估计 : 


P(M|E)P(@) 


PIM) = pop | (7) 
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RJA it, mo 为 次 黑洞 质量 ， q= m2/m1 为 质量 比 ， Xeff 为 有 效 自 旋 参 数 ，Dz 为 光度 距离 。 


图 7 GWTC-2 给 出 了 根据 O3a 观测 到 的 数据 推断 出 的 物理 量 


[11] 
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8 ”更 新 到 O3b 的 GWTC-3 各 个 双 黑 洞 并 合 事件 的 推断 物理 参数 ”的 示意 图 


RF, P(O) 为 参数 的 先 验 分 布 ，P(M) 为 模型 的 证 据 率 ，P(M|9) 为 相应 的 似 然 函数 。 在 
引力 波 数 据 处 理 当 中 ， 可 得 到 ， 


P(M|@) x exp ( CS) et in ») ) (8) 
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其 中 ，s 为 探测 到 的 引力 波 信号 ，n 为 相应 的 高 斯 噪声 ， 圆 括号 代表 信号 相对 于 功率 谱 的 内 
ie 利用 上 述 办 法 ， 则 可 以 对 引力 波源 的 后 验 分 布 做 出 估计 。 由 于 似 然 函 数 非常 复杂 ， 通 

这 个 过 程 需要 利用 马尔 科 夫 链 -蒙特 卡 洛 采样 (MCMC sampling) Be ik Bah (Nested 
e E O 从 而 获取 参数 的 统计 学 范围 估计 。 

为 了 研究 不 同 机 制 产 生 的 双 黑 洞 引力 波 的 自 旋 分 布 ， 可 以 通过 引入 贝 叶 斯 层次 化 模型 
的 方法 来 分 析 集 群 的 信息 。 在 贝 叶 斯 层次 化 模型 中 ， 假 设 不 同 机 制 产生 的 双 黑 洞 引 力 波 事 
件 (FE Nsw) 的 自 施 都 来 自 于 不 同 更 高 层次 对 应 的 集 疼 分 布 ， 而 这 些 分 布 可 以 利用 一 组 参数 
和 进行 描述 区 分 ”。 相 应 的 层次 化 贝 叶 斯 后 验 分 布 可 以 由 以 下 表达 式 表示 : 


p(en Al {h Siig s |Xeg) P D W ， (9) 


RP, iA i 次 观测 到 的 引力 波 信号 (对 应 第 i 个 源 )。 如 果 不 关心 其 中 个 源 的 有 效 自 旋 
分 布 ， 可 以 将 其 边缘 化 ， 得 到 关于 集群 层次 的 参数 分 布 : 


p(Al{a'}) aji Jaap dixis ) plača D) p(A) . (10) 


根据 上 述 公 式 ， 可 以 利用 有 效 自 旋 的 分 布 对 不 同 双 黑 洞 产生 机 制 模型 进行 区 分 和 限制 。 
例如 假设 大 质量 孤立 双星 演化 模型 的 占 比 为 ， 动 力学 机 制 模型 的 占 比 为 1 一 和， 那么 
有 效 自 旋 的 层次 化 分 布 可 以 表示 为 : 


P(xeE|N\) = Pomps(Xer)(1 E 入 ) 十 Paynamical(Xef)A . (11) 


综合 层次 化 贝 叶 斯 方法 和 引力 波 数据 可 以 得 到 关于 模型 占 比 入 的 限制 结 
4.3 ” 自 旋 模型 
前 文 我 们 已 经 提 到 ， a ara ee ed i 
黑洞 的 自 旋 通 常 与 轨道 角 动 量 平 行 ， 即 使 存在 出 生 踢 出 现象 ， 平行 也 不 应 该 被 偏离 得 非 
ie eel Ee ene A 基于 自 旋 分 析 
限制 双 黑 洞 形成 机 制 的 方法 便 是 以 此 为 基础 的 。 通 常 一 般 考 虑 以 下 几 种 不 同 的 模型 ” 。 

(1) DEFAULT 该 模型 假设 双 黑 洞 自 旋 的 大 小 服从 同一 个 Beta 分 布 ， 


T(X1,2 Ax, By) = Beta(ax， By) ) (12) 


其 中 ，ax Alb, A Beta 分 布 的 两 个 参数 。 而 自 旋 倾角 (相对 于 轨道 角 动 量 ) 服从 一 个 包含 两 
种 成 分 泥 合 世 分布: 一 部 分 双 黑 洞 自 旋 与 轨道 角 动 量 准 平行 ( 占 比 为 6)， 对 应 大 质量 双星 演 
化 机 制 的 情况 ， 另 一 部 分 自 旋 则 完全 随机 ， 对 应 动力 学 机 制 的 情况 。 在 平行 的 情况 下 ， 该 模 
型 假设 自 旋 倾 角 服 从 一 个 以 完全 平行 为 期 望 的 截断 正 态 分 布 : 


T (2|¢, 04) = CG: (zlo) +(1—¢)5(z) ， (13) 
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其 中 ，z = cosh 为 双 黑 洞 自 旋 贷 角 的 余弦 值 ，3 为 随机 分 布 ，G (ela) 为 截断 高 斯 分 布 ， 
中 心 为 > = 0，o 为 标准 差 。 我 们 可 以 看 到 这 样 的 模型 从 直观 上 最 简洁 ， 一 共 只 包含 4 个 参 
数 (ax, Bx, Ç, O1)e 

(2) GAUSSIAN 该 模型 的 优点 在 于 可 以 用 一 个 联合 高 斯 分 布 同时 拟 合 Xer 和 X 两 


DHE: 
T (Xer, Xpl| Hef, Oeff, Hp, Op, p) xG (Xef， Xp|H， 5) ， (14) 
其 中 ， H= (Hef, Hpo 协 方差 矩阵 为 : 
O25 Daefa 
y= $ 5 p ; (15) 
POeff Op os 


该 模型 总 共有 5 个 参数 (Metts Com, Hp: Tp, P) MM oeg Mop 将 会 帮助 我 们 分 析 自 旋 诱 导 的 轨 
道 进 动 效应 以 及 部 分 具有 反 平 行 的 自 旋 的 双 黑 洞 系 统 。 

(3) MULTI SPIN a D a 布 和 目 旋 分 布 ， 其 参数 众多 (12 
个 自 旋 参 数 ，10 个 质量 参数 )， 主 要 用 于 分 析 自 旋 - 质 量 、 质 量 比 关 系 。 这 种 模型 通常 被 用 于 
分 析 大 质量 的 双 黑 洞 系 统 是 否 eee 于 等 级 多 体 的 并 合 。 

4.4 当前 结果 

本 节 ， 我 们 将 简要 总 结 基于 LIGO/Virgo 引力 波 探测 获得 的 双 黑 洞 自 旋 分 布 限于 

洞 产 生机 制 的 相关 结果 和 进展 。 

4.4.1 GWTC-1 及 以 前 

自 GW150914” 被 发 现 ， 便 立即 有 一 些 学 者 对 其 形成 机 制 进行 了 研究 。 例 如 ， 一 些 学 者 
认为 40Mo ~ 100Mo 区 间 的 双 黑 洞 引力 波源 可 以 由 孤立 双星 演化 过 程 形成 ”， 同 时 也 有 一 
些 学 者 认为 动力 学 演化 也 可 以 形成 此 次 观测 到 的 双 黑 洞 “。 在 GWTC-1 前 虽然 观测 到 的 
引力 波源 并 不 多 ， 人 们 已 经 开始 开展 波源 竹 质 和 相关 模型 的 统计 学 研究 ， 特别 是 利用 自 旋 分 
布 分 析 双 黑洞 形成 机 制 的 探究 。 例 如 ， 通 过 对 有 效 自 旋 xe 分布 简单 贝 叶 斯 分 析 ” 发现， 如 
果 认 为 探测 到 的 黑洞 的 自 旋 较 大 ， 那 么 自 旋 方向 的 分 布 更 加 倾向 于 各 向 同性 ; 而 在 小 自 族 的 
先 验 条 件 下 ， 这 一 倾向 性 将 会 大 大 减弱 。 图 加 给 出 了 各 向 同性 分 布 在 双 黑 洞 自 旋 后 验 分 布 
中 可 能 的 占 比 。 值 得 注意 的 是 ， 由 于 当时 的 观测 数据 很 少 ， 所 以 无 法 在 统计 上 说 明 这 种 结论 
是 普遍 的 。 

44.2 GWTC-2 

利用 上 述 的 DEFAULT 模型 ， 发 现 了 GWTC-2 siege 自 旋 -轨道 角 动 量 的 非 
平行 现象 ， 见 图 [LU 这 一 现象 的 发 现 意味 着 一 部 分 双 黑 洞 要 么 具有 完全 随机 的 自 旋 方向 ， 
要 么 黑洞 的 自 旋 倾 角 几 乎 为 0 但 是 具有 一 定 的 弥散 。 ala 种 图 景 ， 都 表明 了 某 些 双 
黑洞 存在 轨道 角 动 量 平面 内 的 geld 图 10 展示 了 利用 DEFAULT 模型 得 到 的 双 黑 洞 自 
旋 大 小 以 及 自 旋 倾角 的 后 验 分 布 ， 其 = |s| 为 无 量 纲 自 旋 大 小 。 

虽然 GWTC-2 Meee eer a 果 ， 但 是 其 中 的 弥散 也 不 能 忽略 。 此 
外 ， 如 果 利 用 GAUSSIAN 模型 ， 可 以 在 99% 的 置信 度 上 排除 up = op = 0 的 情况 ， 这 表明 


= 


双 黑 


= 
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一 后 验 Flat 
-'- 后 验 Low 


0.0 0.2 0.4 0.6 0.8 1.0 
各 向 同性 占 比 £ 


VE: 其 中 的 High, Flat. Low 分 别 表示 假设 无 量 纲 自 旋 大 小 的 分 布 分 别 取 2ay. 1, 2(1 一 as) 时 孤立 双星 演 
化 的 后 验 分 布 。a 代表 无 量 纲 自 旋 大 小 范围 为 0 到 1。 


9 ”GWTC-1 对 动力 学 以 及 大 质量 恒星 演化 模型 占 比 的 预测 分 布 ” 


x 
a) 


VE: 图 中 的 橙色 亮 线 代 表 每 个 独立 观测 事件 下 给 出 的 无 量 纲 自 旋 幅 度 的 分 布 ， 黑 线 给 出 无 量 纲 自 旋 幅 度 的 分 布 ， 

虚线 代表 90% 置信 度 区 域 。 
图 10 a) 给 出 了 由 DEFAULT 模型 推断 下 的 无 量 纲 自 旋 幅度 x 的 分 布 关系 ; b) 给 出 的 是 重 构 后 两 个 黑 
洞 自 旋 与 双 黑洞 轨道 角 动 量 平面 夹 角 cos 9, 的 分 布 关系 


自 旋 诱导 的 轨道 进 动 现象 的 存在 ， 也 从 另外 一 个 侧面 验证 了 自 旋 与 轨道 角 动 量 的 非 平行 现 
象 的 存在 。 

其 次 ， 人 们 在 GWTC-2 的 双 黑 洞 样本 中 发 现 了 具有 反 平 行 的 例子 。 这 提示 可 能 有 不 止 
一 种 双 黑 洞 的 产生 机 制 。 利 用 GAUSSIAN 模型 推测 出 某 些 系统 可 能 存在 Xer < 0 的 现象 ， 
这 等 价 于 双 黑 洞 中 至 少 有 一 个 与 轨道 角 动 量 的 夹 角 大 于 90°. 

图 四 展示 了 边缘 化 op 和 Xp 之 后 平均 per 和 oer 的 后 验 分 布 ， 峰 值 分 别 在 0.06+008 
和 0.12+0:04。 这 要 求 部 分 双 黑 洞 的 有 效 自 旋 Xer < 0。 而 利用 DEFAULT 模型 ， 我 们 也 可 
以 得 到 相似 的 结论 。 图 11 a) 图 展示 了 利用 该 模型 重建 的 cosh 的 分 布 ， 其 中 可 以 看 到 
cos 91,2 < 0 的 可 能 性 仍然 存在 ， 这 意味 着 反 平 行 的 自 旋 现象 的 存在 。 
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0.06 


11 ”利用 GAUSSIAN 模型 给 出 双 黑 洞 等 效 自 旋 ver 的 期 望 值 es 与 标准 差 oog 的 分 布 


利用 图 四 展示 的 er 和 oon 的 后 验 分 布 ， 可 以 获得 关于 户 (Xer > 0)、f (Xer <0)» fo 
(ver = 0) 的 占 比 后 验 分 布 ( 见 图 ED。 可 以 看 到 在 99% 的 置信 和 度 的 情况 下 ，f, = 0.672228, 
fa = 0.274917 Ful fy = 0.05+0.02。 就 像 之 前 提 到 过 的 ， 动 力学 机 制 并 不 是 产生 某 些 双 黑 洞 具 
小 于 零 的 有 效 自 旋 现象 的 唯一 方法 。 如 果 前 身 星 在 超新星 爆发 阶段 经 历 了 一 个 很 强 的 出 
E 踢 出 和 几乎 无 效 的 自 旋 再 平行 过 程 ， 由 大 质量 孤立 双星 演化 机 制 产生 的 双 黑 洞 中 约 小 于 
0% 可 能 具有 小 于 0 的 有 效 自 旋 。 这 之 中 还 存在 很 多 复杂 性 ， 例 如 在 孤立 等 级 三 体系 统 中 ， 
岂可 能 产生 具有 非 平行 自 旋 分 量 的 双 黑洞 。 

尽管 如 此 ， 如 果 简单 地 假设 xer < 0 的 双 黑 洞 均 产 生 于 星团 中 的 动力 学 相互 作用 ， 则 可 
以 利用 fa 来 限制 动力 学 起 源 的 双 黑 洞 的 比率 ， 那 么 相应 的 由 动力 学 机 制 产生 的 双 黑 洞 的 占 
比 为 : 


Ee it 


fn 
a ae (16) 
以 及 孤立 双星 机 制 的 占 比 为 有 


在 90% WEARER, fa 的 值 在 0.25~0.93 之 间 。 这 说 明 GWTC-2 中 的 双 黑 洞 一 部 分 通过 
抓 立 双星 机 制 形 成 ， 另 一 部 分 则 是 通过 星团 中 的 动力 学 机 制 产 生 。 关 于 动力 学 与 孤立 双星 演 
化 模型 下 讨论 有 效 自 旋 参 数 的 产生 细节 以 及 为 何 采用 上 述 参数 描述 有 更 详细 的 讨论 “ ”。 

此 外 ， 如 果 GWTC-2 数据 中 存在 等 级 并 合 的 现象 ， 自 旋 将 会 与 双 黑 洞 系统 的 质量 存在 
一 定 的 关系 ， 即 越 重 的 等 级 双 黑 洞 并 合 将 会 拥有 更 大 的 自 旋 。 但 是 利用 上 文 提 到 的 MULTI 
SPIN 模型 对 GWTC-2 数据 进行 分 析 ， 并 没有 发 现 自 旋 分 布 与 质量 存在 关联 的 证 据 。 这 一 
点 与 GWTC-1 得 到 的 结果 是 类 似 的 。 


0.277015 


0.05701 


12 ”利用 GWTC-2 BGR S3IM fo. fo fo MEMO” 


44.3 GWTC-3 

KALE GWTC-2, GWTC-3 双 黑 洞 样本 扩大 了 很 多 ， 但 是 利用 这 些 双 黑 洞 的 自 旋 分 布 
统计 分 析 得 到 的 很 多 结论 都 非常 相似 ， 如 图 苔 展示 了 利用 DEFAULT 模型 对 双 黑 洞 的 自 旋 
大 小 以 及 自 旋 倾角 的 限制 。 可 以 看 到 ， 相 比 GWTC-2，GWTC-3 双 黑 洞 自 旋 倾角 的 余弦 值 
cos9 分 布 更 加 平坦 ， 更 加 支持 具有 完全 随机 的 自 旋 倾角 分 布 。 又 如 ， 利 用 GWTC-3 双 黑 洞 
样本 ， 由 GAUSSIAN 模型 可 以 得 到 xer 和 x, 的 后 验 分 布 ， 发 现 wee JES ( 见 图 14); 这 一 
结论 与 GWTC-2 的 结果 一 致 。 

根据 GWTC-3 的 数据 得 出 了 两 个 非常 重要 的 新 结论 。 一 是 自 旋 的 分 布 在 不 同 的 质量 范 
围 都 是 一 臻 的。 图 四 展示 了 自 旋 与 轨道 角 动 量 的 平行 分 量 s。 的 大 小 和 双 黑 洞 系统 的 咽 咯 
质量 的 关系 。 在 小 质量 端 ， 平 行 分 量 几 乎 为 0; 而 在 大 质量 段 ， 即 使 存在 比较 大 的 弥散 ， 平 
行 分 量 的 大 小 分 布 仍然 和 0 一 致 。 对 咽 嗽 质量 小 于 30M。 的 双 黑 洞 自 旋 平行 分 量 取 平 均 为 
0.38 (90% 置信 度 )， 而 在 小 质量 端 该 结果 为 0.5。 

另外 一 个 是 在 GWTC-3 的 数据 中 发 现 了 质量 比 qg 与 有 效 自 旋 Xeg 的 反 相 关 关 系 ， 因 此 
质量 比 趋 近 于 1 的 双 黑 洞 系 统 倾向 于 拥有 接近 0 的 有 效 自 旋 xeg， 而 质量 比 很 大 的 系统 则 倾 
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---- GWTC-2 
oe Cee oe —GWTC-3 


-0.6 02 00 02 
Xett Xp 
a) b) 
图 14 ”数据 更 新 到 CWTC-3 后 有 效 自 旋 xe URS xo HERSE” 
30 40 50 60 
M/M: 
图 15 


在 GWTC-3 现 有 数据 下 给 出 的 |s-| 与 咽 嗽 质量 的 限制 关系 
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FF 
布 p( Xe): 


其 中 ， 


洞 的 有 效 自 旋 大 于 0 的 可 


图 16 


5 


引力 波 的 发 现 极 大 地 拓展 了 人 类 对 宇宙 


(gq) 和 o(a) 用 于 


I 有 一 个 较 大 的 正 有 效 自 旋 。 利 月 


K Pe - 


ERK. B 


天 文学 


(Xer = 
20(q) 

HIB Xef 分 布 的 中 值 和 标准 差 随 q 的 变化 关系 。 

叶 斯 方法 ， 可 以 在 97.5% 置信 度 下 认为 ula) 为 9 的 减 函数 ， 这 也 意味 着 质量 比 越 大 ， 双 黑 

ERRI q 与 Xeg 的 后 验 分 


进 展 


HL(d)) 


E 


41% 


H GAUSSIAN 模型 ， 可 以 获得 Xor XT q 的 条 件 后 验 分 


(18) 


利用 上 述 层 次 化 贝 


) 布 。 


1.0 


0.6F 


总 结 吉 与 展望 


0.00 
x 


eff 


PC AER A AY 


0.50 


有 效 自 旋 参 数 与 质量 比 在 GWTC-3 的 观测 结果 中 呈现 出 的 反 相 关 的 关系 示意 图 


E 解 与 认识 。 


1.00 


[41] 


引力 波 探 测 到 


通过 


大 量 恒星 级 双 黑 洞 并 合 ， 其 究竟 是 如 何 形 成 并 演化 至 最 后 的 并 合 仍 是 未 解 之 谜 。 人 们 提出 了 


多 种 形成 双 黑 洞 的 机 制 ， 


这 些 过 程 的 理解 还 


精准 


不 同 机 种 


本 中 自 旋 -轨道 角 动 量 的 非 平和 
的 现象 等 。 这 对 限 


化 机 制 以 及 倾 


随 着 引力 波 探测 的 持续 
的 第 三 代 地 基 引 力 波 探测 器 ， 例 如 爱 
(Cosmic Explorer, CE) 将 于 2030 FARZ 


Rel] XX 


所 其 起 源 ， 尤 其 是 
向 于 各 向 同性 随机 分 布 的 动力 学 本 
开展 ， 人 们 将 更 准确 


区 分 双 黑 洞 


自 旋 等 物理 参数 的 分 布 上 有 明显 的 区 别 。LIGO/Virgo $ 
的 测量 为 人 们 提供 了 一 个 PB 


HHEH HL 


上 在 双 黑 洞 形成 进程 中 都 包含 了 众多 天 体 物理 过 程 ， 
很 大 的 不 确定 性 。 尽 管 如 此 ， 人 们 普遍 相信 不 同 机 制 产 生 的 双 黑 洞 在 质 
等 地 基 引 力 波 探测 器 对 自 旋 较为 
洞 不 同形 成 机 制 的 绝 佳 手段 。 
斯 方法 ， 利 用 现 有 的 GWTC 数据 集 ， 人 们 已 经 取得 一 些 初步 的 进展 。 


目前 对 


基于 自 旋 分 析 的 贝 叶 
网 如 发 现 了 双 黑 洞 样 


了 现象 ， 自 旋 与 双 黑 洞 质量 几乎 无 关 ， 而 与 其 质量 比 呈 反 相 关 
两 成 分 自 旋 方向 倾向 于 平行 的 大 质量 双星 演 
Hl, BATA 


深刻 的 意义 。 


eig 


Sh Ml) ee Te Fe BQ MY A 


& 洞 自 旋 分 布 。 特 别 是 未 来 


因 斯 坦 望 远 镜 (Einstein Telescope, ET) 和 宇宙 勘探 者 


至 渐 建成 。 它 们 极 高 的 探测 灵敏 度 将 会 大 大 地 


扩展 恒 


4 期 陆地 ， 等 :恒星 级 双 黑 洞 自 旋 研究 进展 473 


星 级 双 黑 洞 的 引力 波 探测 样本 和 改善 物理 参数 测量 的 精度 。 随 着 数据 的 不 断 积累 ， 关 于 双 黑 
洞 的 自 旋 统 计 学 分 析 将 会 帮助 人 们 最 终 揭 示 恒 星 级 双 黑 洞 的 起 源 ， 有 望 限制 黑洞 形成 过 程 
中 的 众多 重要 物理 过 程 并 揭秘 黑洞 自 旋 的 产生 。 


H 审 稿 人 对 此 文章 提出 的 详细 修改 意见 ， 同 时 感谢 曾 思 藻 对 本 文 攻 
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Progress in Constraining the Origin of Stellar Binary Black 
Holes via Their Spin Distribution 


LU Di!?, CHEN Zhi-wei!:?, YAN Chang-shuo!?, LU You-jun!? 


(1. National Astronomical Observatory, Chinese Academy of Sciences, Beijing 100101, China; 2. School 
of Astronomy and Space Sciences, University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: Since the first detection of GW150914, a larger amount of stellar binary black 
hole (sBBH) mergers have been detected by LIGO/Virgo collaboration with high confidence. 
By analyzing the gravitational wave (GW) signals, one can obtain crucial physical properties 
of these sBBH mergers, including their chirp mass, spin, and luminosity distance, etc. The 
origin of sBBHs may encode in the distribution of these physical properties. Here, the 
curent status of the study on sBBHs spins is reviewed. First, various formation mechanisms 
of sBBHs are briefly introduced, and these mechanisms include: 1) sBBHs are formed via the 
evolution of massive binary stars (EMBS channel); 2) sBBHs are formed via the dynamical 
interactions with surrounding (compact) stars in dense stellar clusters (dynamical channel); 
etc. The origin of spins for sBBHs originated from the black hole formation processes is 
also introduced. One may expect that the spin distributions resulting from different sBBH 
formation channels are different. For example, the two components of each sBBH resulting 
from the EMBS channel have spins more or less aligned with each other, while they may 
be randomly distributed if the sBBH is resulting from the dynamical channel. Therefore, 
the spin distribution obtained from GW observations may be used to distinguish the EMBS 
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channel from the dynamical channel and further reveal the origin of sBBHs. The general 
methods for constraining the origin of sBBHs via their spin distribution are then introduced. 
The results obtained by recent works using the GW observations of sBBHs (GWTC) are 
summarized. Finally, a summary and a short prospect are given for the future demography 


study of the sBBH spins measured by next generation ground-based GW detectors. 


Key words: gravitational wave astronomy; gravitational wave Sources; black holes 


